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Abstract 

We propose a cosmological model which could explain, in a very natural way, 
the apparently periodic structures of the universe, as revealed in a series of recent 
observations. Our point of view is to reduce the cosmological Friedman-Einstein 
dynamical system to a sort of Schrodinger equation whose bound eigensolutions 
are oscillating functions. Taking into account the cosmological expansion, the large 
scale periodic structure could be easily recovered considering the amplitudes and 
the correlation lengths of the galaxy clusters. 
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1 Introduction 


Understanding the large scale structure of the universe is one of the main issue of modern 
cosmology. In fact, a large amount of observations tell us that the cosmological princi¬ 
ple breaks at scales of the order ~ 100/i“^Mpc and below, while, well above this limit, 
(~ 1000/i“^Mpc) it seems that the matter in the universe is homogeneously and isotrop¬ 
ically distributed (Kolb & Turner 1990). The symbol h indicates the normalized Hubble 
constant whose values are 0.4 < h < 1. 

From a cosmological point of view {i.e. we are considering scales where galaxies are 
pointlike constituents), the clustering properties of objects clearly show that a sort of 
hierarchy exists: we have tight and loose groups of galaxies, clusters, superclusters and 
filaments of galaxies. All these structures come from an excess of correlation over given 
distributions and they can be well represented by the two-point correlation function 

e,(r) = , (1) 

which hxes the correlation scales (Kolb & Turner 1990). The index j can be j = G,C, S 
if the correlation is between galaxies, clusters, or superclusters. Cj determines the cor¬ 
relation amplitude which can be from 0.1 to 10/i“^Mpc for galaxies, up to 25/i“^Mpc 
for clusters, and up to 100/i“^Mpc for superclusters (Peebles 1993). The exponent pj 
assignes the slope of the power law and it is found that, for all the hierarchical orders, 
it is the same, that is nj ~ 1.8. In other words, the excess of correlation is the “same” 
for galaxies, clusters and superclusters. This fact could mean that the distribution of 
objects is absolutely non-random and that some fundamental mechanism has led to the 
formation of large-scale structure. 

A further evidence of such a point of view is the result by Broadhurst et al. (from 
now on BEKS) (Broadhurst et a/. 1990) which found an excess of correlation and an 
apparent regularity in the galaxy distribution with a characteristic scale of 128/i“^Mpc. 
The structure was revealed after the completion of deep pencil beam surveys extending to 
red-shift 2 : > 0.2 toward northern and southern Galactic poles (Broadhurst et a/. 1990). 
The linear extension of the surveys was to 2000/i“^Mpc so that, the scale of 128/i“^Mpc 
between successive density peaks could be interpreted as a sort of transition scale to 
homogeneity and isotropy at larger scales. 

Several interpretations have been given to this seminal result. For example, cosmo¬ 
logical solutions with an overall Friedmanian expanding behaviour, corrected by a small 
oscillatory regime, could implement such a structure (Busarello et a/. 1994). Furthermore, 
dynamics of a nonconformal scalar held could modify cosmological evolution in the sense 
of oscillations (Morikawa 1991), (Capozziello et a/. 1996). 

It emerges that several observational quantities can be affected by such a behaviour 
since oscillations in the cosmological solutions can be easily reduced to oscillations in the 
red-shift ; 2 . In fact, if a{t) is the cosmological scale factor (function of the cosmic time 
t), we have the relation 



so that all quantities containing the Hubble parameter H or the red-shift have to 
oscillate. Then the oscillation in red-shift can be considered, as first proposed by Tifft 
(Tifft 1977, Tifft 1987), as a sort of “quantization”. On the other hand, the apparent 
periodicity in the distribution of galaxies implies peaks in the red-shift distribution which 
lay on concentric spherical shells with periodically spaced radii (someone claimed for a 
new Tolemaic principle from the point of view of Earth!). 

The main quantities affected by oscillations are the following. 


1. The number count-red-shift relation of galaxies 


dN 

dVLdLdz 


n{L,to)alH , 


( 3 ) 


where dN is the number of galaxies in the solid angle dO having red-shift between 
and z + dz and luminosity between L and L + dL] n{L, to) is the number density 
of galaxies with luminosity L that an observer sees at to] Uq is the actual scale 

factor and d is the comoving distance dehned as d = / —. If id oscillates, the 

Jt a 

number of galaxies in a solid angle changes as a function of the red-shift and the 
luminosity. 


2. The two-point correlation function which can be rewritten in term of H as 


= 



( 4 ) 


where r is the red-shift distance (Weinberg 1972) 


r = zHq ^ 



( 5 ) 


with L the luminosity of the object and F the measured flux; ds is the comoving 
distance of the sample, is the sample size, which can be r* ~ 5h“^Mpc for the 
galaxies and ~ 25/i“^Mpc for the clusters; Ho is the actual Hubble constant. 
However, the oscillations can enhance or depress the galaxy-galaxy correlation 
function since the distance between the objects apparently increases or decreases. 


3. The number count-redshift relation for quasars. Since the oscillations of the Hubble 
parameter is coherent and global, we expect a similar clustering property also in 
higher red-shift regions such as 1 < < 6. 


In general, given a sample of homologous objects spaced in red-shift, we can, in princi¬ 
ple, apply the same argument used for the galaxies, that is we can expect an apparent 
periodicity in the distributions [e.g. Lyman a-clouds, radio galaxies, and so on). 

Besides we can ask for some fundamental physics explaination of such a phenomenol¬ 
ogy as the above mentioned dynamics of a scalar held. For example, some authors are 
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recently supposing a deep link between quantization and gravitation. Some of them 
support the point of view that gravity could give rise to a sort of “stochastic” quantiza¬ 
tion which relates the cosmic scales of interest to the Planck constant (Calogero 1997). 
Others search for extensions of Einstein’s relativity able to include, for examples, fractal 
structures (Nottale 1997). A more conservative point of view could be considering cosmic 
structures as ’’remnants” of primordial quantum processes and trying to connect early 
“quantum” dynamics with today-observed classical macroscopic dynamics. 

This is, in some sense, the usual approach of most of quantum gravity theories which 
consider either the whole universe as a quantum system (supposing, for example, several 
co-existing, noninteracting universes (Everett 1957, De Witt 1967, Vilenkin 1982) or the 
universe as a classical background where primordial quantum processes have given rise 
to the actual macroscopic structures (this is the approach of quantum held theory on 
curved spacetimes (Birrell & Davies 1982). 

In relation to the hrst point of view, it has been shown that the Wheeler-DeWitt 
equation 'H'lp = 0 for the wave function of the universe ip can be written, in the 

dip 

Friedman-Robertson-Walker (FRW) spacetime, as a Schrodinger equation = Hip, 

da 


where p, 


da 

dt 


(Pollock 1997). In this sense {i.e. in the mini-superspace approximation). 


quantum gravity can be treated as ordinary quantum mechanics. This approximation 
could imply experimental verihcations of the Wheeler-DeWitt equation (Pollock 1997, 
Colella et a/. 1975) and give further constraints on big-bang nucleosynthesis (Llorente & 
Perez-Mercader 1995). 

In this paper, following Rosen 1993, we want to reduce, at least formally, the cos¬ 
mological Einstein-Friedman equations of general relativity to a quantum mechanical 
system. If this issue holds, the Friedman equations can be recast as a Schrodinger equa¬ 
tion and the cosmological solutions can be read as eigensolutions of such a “cosmological” 
Schrodinger equation. In this quite simple scheme, the bound states (with negative en¬ 
ergy) are automatically oscillating functions and the probability amplitudes can £t the 
BEKS oscillatory distribution with an opportune choice of the parameters: that is, it is 
nothing else but a particular eigensolution. This means that the observed oscillations are 
just the remnant of a primordial quantum state enlarged by the cosmological expansion 
(which has been huge if an early inflationary phase is supposed). An important point 
has to be stressed: here, we are not considering the above mentioned minisuperspace ap¬ 
proach in which one estimates the probability that a particular “classical” universe comes 
out from quantum boundary conditions (Everett 1957, De Witt 1967, Vilenkin 1982) : 
we are just recasting the cosmological dynamical system as a Schrodinger quantum-like 
system where bounded solutions could give rise to the observed large scale structure. The 
probability has to be connected to the above correlation function. 

The paper is organized as follows. In Sec.2, following the approach by Rosen (Rosen 
1993), we construct our cosmological Schrodinger equation. Sec.3 is devoted to solve such 
an equation stressing the solutions characterized by oscillatory behaviours. Discussion of 
the results and conclusions are drawn in Sec.4. 
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2 Cosmological equations as a Schrodinger equation 

The cosmological evolution equations come from the Einstein held equations 


4 Rfiu T ) (6) 

^ C 

where is the Ricci tensor, R is the Ricci scalar, is the metric tensor, is the 
stress-energy tensor, and A is the cosmological constant which is taken into consideration 
in order to generalize the discussion. Such equations can be derived from the Hilbert- 
Einstein action 

kl = — \QqYQ j [R + 2 A + Crn] , (7) 

where Cm is the Lagrangian density relative to the matter helds. As a simple hypothesis, 
we consider standard perfect huid matter so that, by varying Cm with respect to we 
get 

= (p + p)u^Uy - pg^u , (8) 

where p and p are the energy-matter density and the pressure respectively. The dynam¬ 
ical problem is completely set when matter evolution equations are given; they are the 
contracted Bianchi identities 

= 0 . ( 9 ) 

A further equation has to be imposed in order to assign the thermodynamical state of 
matter. It, usually, is 

P = (7 - l)cV , (10) 

where 7 is a constant (1 < 7 < 2 for standard perfect fluid matter). 

The Einstein-Friedman cosmological equations are easily recovered as soon as the 
Friedman-Robertson-Walker metric 

ds^ = c^dt^ — a^{t) I ——7VT + r'^dO'^ + sin^ 9d(i)^\ , (11) 

(1 — kr^ J 

is imposed. Here a(t) is the scalar factor of the universe, k = 0, ±1 is the curvature 
constant. 

Then, from Eqs.(|^) and (P), we obtain the system 


d 47rG , , Ac^ 

—- 3 ^(^> + 3rt+ 3 . 

(12) 

/ d\^ kc^ 8nG Ac^ 

\a/ 3c^ ^ ' 3 ’ 

(13) 

P + 3 (p + p) = 0, 

(14) 


which has to be completed by Eq. (pT|) . 
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( 15 ) 


Let us now rewrite the (0, 0)-Einstein equation (|^ as 


-ma^ 


m /SttG 


m / 


6 V 


p + , 


where m is a mass which will be specihed below. 

We can formally write, following Rosen (Rosen 1993) 


T + V = E, 


(16) 


where 


and 


T = 


V{a) = 


1 
2 

dvrmG 


3c2 


pa? — 


mKc^a? 



(17) 

(18) 

(19) 


T has the role of the kinetic energy, V{a) of the potential energy, and E of the total 
energy of the system. 

Since Eq.(|^) is the hrst integral of Eq.(^), the equation of motion of a “particle” of 
mass m is recovered, being 

dV , , 

ma = —— , (20) 

da 

and using Eq.(|^. 

By this formalism, we can dehne a “particle” momentum 


n = md , 


and the Hamiltonian 




Using the simple scheme of hrst quantization, we can substitute 

. d 

n , 

da 

and then derive the Schrodinger equation 


where the wave function is 


( 21 ) 

( 22 ) 

( 23 ) 

( 24 ) 

( 25 ) 
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In our case, m can be interpreted as the mass of a galaxy and Itbp as the probability to 
hnd such a galaxy at a given a{t). Immediately, we can pass from a{t) to the observable 
red-shift z by the relation 

— = l + z, (26) 

a 

where oq is the actual scale factor. In this sense, tb = is the formal probability 

amplitude to hnd a given object of mass m at a given red-shift z, at time t. 

Following the ordinary quantum mechanics, a stationary state of energy E is 

^(a,f) = (27) 

and the Schrodinger stationary equation is 

-^p^ + V(,a)i, = Ei,, ( 28 ) 

2m aa^ 

or, in terms of z, 


r (z + 1 )^ (P'4) 2{z + 1 )^ djj 

2m Oq dz‘^ Og dz 


+ V{z)tl: = E'lP , 


(29) 


being 


V{z) 


4:7rmG\ p{z)al mAc^ 
3 c ^ / (1 + ^)^ 6(1 + z )2 


(30) 


We shall continue our analysis by using the scale factor a{t), since the results can be 
easily translated in 2 ;. 

To conclude this section, we have to stress that our model differs from that in Rosen 
( Rosen 1993), since there m was the mass of the whole universe, while we are using an 
effective “particle”, which, from a cosmological point of view, can be a galaxy. Further¬ 
more, we propose a different probability interpretation for T. 


3 The eigensolutions 

From Eqs. (0 and (p!^, we obtain 


— An-^z 


p = Aa 


( 31 ) 


where R is a positive integration constant which gives the density of standard matter 
at a given epoch og. The thermodynamical state of matter is assigned by giving 7 . We 
have, for example, 7 = 1 for dust-like matter and 7 = 4/3 for radiation. 

The Schrodinger stationary equation (^) can be written in the form 






V’ 


0 , 


(32) 
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where 


8nGm^ rrrc^A 2mE mc^k 

and the prime indicates the derivative with respect to a. Snch a parametric differential 
eqnation can be exactly solved in several cases giving rise to oscillatory behaviours. 
In any case, being nothing else but the one-dimensional problem of ordinary quantum 
mechanics, we can discuss its general features, in particular when oscillatory solutions 
are present. 


3.1 General considerations 


Let us rewrite Eq.(^) as 


where 


v-" +[S- 17(a)] = 0 , 


U(a) 




(34) 

(35) 


By comparing the potential 17(a) and the “energy” £, we can discuss the qualitative 
behaviour of solutions. However, the form of potential strictly depends on the parameters 
H, 7 , and C (alternatively on A, 7 , and A) so that, due to the meaning which we are 
going to give to ^fJ, the clustering and correlation properties are functions, as it must be, 
of the kind and of the content of the fluids involved into the dynamics. Furthermore, the 
energy is ’’fixed” by the given FRW-spatial model. 

If we require that we cannot find galaxies at infinity, the function 'ijj{a) is normalized. 
As a trivial result, the function 'ijj{a) = 0 is a solution of the problem. To avoid that this 
is the only solution, we must ask for nontrivial normalized solutions which will set up a 
discrete spectrum. 

By fixing the energy (that is by choosing the FRW model), we can divide the a-axis 
in a certain number of intervals using the points where T — 17(a) =0. We assume, for 
simplicity, that such points are a finite number. Then, the extreme intervals have an 
unbounded extension, and, between them, we have a finite number of bounded intervals. 
We call type-/ intervals those where T < 17(a) so that '0'VV’ > 0; and type-// intervals 
those where S > U{a), so that iIj"/ il) < 0. In type-/ intervals, the function -0 never 
changes its sign, in type-// intervals, it is, in general, oscillating and can change it sign. 
In this case we could have solutions showing clustering and anti-clustering properties for 
galaxies. We have bending points if T = 17(a): in these cases, the solution transits from 
type-/ to type-// intervals and viceversa. The function is normalizable if the extreme 
intervals are of type I. In fact, if we are considering type-// intervals, the “particle” 
can escape towards infinity, that is it can stay at extremely great distance from a given 
point. This fact is physically in contrast with any clustering or correlation property. In 
other words, the discrete spectrum is inside the energy values which do not allow the 
particle to escape towards infinity. For this reason, the energy must be less than 17(a). 
However, we have to stress that the situation t —>■ —00 means that we are considering 
singularity-free models but a{t) has to remain nonnegative. 
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If we are in some regime where \£\ S> U{a), with £ positive dehned (z.e. /c = — 1), 
Eq. (^H) has stationary wave solutions of the form 

V'(a) = + 026 "*^“ , (36) 

if the energy is negative, the physical solutions are decreasing exponential functions. Such 
functions are not always normalizable but it is easy to see that oscillatory behaviours 
can be recovered. 


3.2 Special cases 

From Eq.(^4|), by choosing suitable values of the parameters S, 7 , C, and £, we can 
recover all the standard cases of ordinary quantum mechanichs (potential well, harmonic 
oscillator. Coulomb-like potential, and so on). As it is well known (Messiah 1961), many 
of them exhibit oscillatory eigensolutions. Considering our point of view, the interesting 
values of 7 are 7 = 4/3 (radiation), 7 = 1 (dust), 7 = 0 (scalar held). 


3.2.1 The radiation case 

In this case, Eq. (|3^) becomes of the form 

B 


r + 


+ Ca^ + £ 




-0 = 0 


For A = 0 and /c = — 1, it can be written as 

,2 - 1/4) 


r + 






^ = 0 , 


with 


(3^ = £ = 


me 


2 ’ 


v = \ - B 


(37) 


(38) 


(39) 


which is a Bessel equation. The general solution is (Abramowitz & Stegun 1970) 

-(/(a) = y/aZ^ {/3a) . (40) 


From the theory of Bessel functions, Zy{x) can be given as a combination of the func¬ 
tions Jy{x), Yy{x), or H‘^\x). The analysis of asymptotic behaviours results 

particularly interesting. We have, for a —0, 


jU ((^) 
Y,{a) 

YM 


(i)T(..,, 


-iH^\a) ^ iH'^^a) ~ — Ina, 

TT 

-iHl^\a) ~ iH’3){a) ~ -^r(z/) 


r(2)/ 


r\j 


(41) 









In the opposite situation, that is a —cx), we have oscillating solutions 


J i> (^) 

YM 

Hf\a) 

For A = 0 and k = +1, we get 



( 42 ) 

( 43 ) 



( 44 ) 


where the relations 


Jy{ia) = Iy{a ), 

Y,{ia) = , 

TT 

H^^\ia) = (45) 

TT 

hold. 

If /c = 0 but A 7 ^ 0, we have 

i){a) = ^/aZ^ 


where 


VI-45 


( 47 ) 


3.2.2 The dust case 

In the dust-dominated case, we have 


v" + 


- + Ca‘^ + S 
a 


V = 0, 


and for k = A = 0, we have solutions of the form 

'ipia) = A/aZ.,(\/a), 


( 48 ) 


( 49 ) 


where the above considerations hold. If C* = 0, we get an hydrogenlike Schrodinger 
equation (with null angular momentum) and the wavefunction V(®) can be written as a 
combination of Laguerre polynomials (see also Rosen 1993, Messiah 1961). 
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3 . 2.3 The scalar field case 


The scalar field case, is not exactly a standard perfect fluid form of matter but it is 
of extreme physical interest since it is connected to inflation (Kolb & Turner 1990). 
However, it can be treated as a standard fluid assuming 7 = 0 in the equation of state 
Eq.([To|) (Guth 1981). From the point of view of dynamics, it has the same role of 
cosmological constant so that Eq.(^) becomes 




Ca^ + £ 


-0 = 0 


( 50 ) 


We get the same equation by putting B = 0. However, the cosmological constant must 
be redefined. It is clear that, depending on the relative signs and values of C and S, 
Eq.(^^ reduces to the equation of an harmonic oscillator with combinations of Hermite 
polynomials as eigensolutions or, in general, to the differential equation of parabolic 
cylinder functions (Abramowitz & Stegun 1970). 


4 Discussion and Conclusions 

Several theoretical interpretations of BEKS’s result have been attempted so far (Busarello 
et a/.1994, Morikawa 1991, Capozziello 1996, Budinich & Raczka 1993) but in all of them 
modifications of Einstein standard theory of gravity have been invoked. Essentially, a 
further ingredient as a scalar field should have had the role of perturbing the overall 
Friedmanian dynamics of the universe breaking the homogeneity and isotropy at scales 
of the order of ~ 100h“^Mpc. The approach which we have adopted in this paper is, 
in some sense, absolutely conservative since we are still using Einstein general relativity. 
The only issue is that the Friedman-Einstein cosmological dynamics could be read in 
the sense of elementary quantum mechanics: the ( 0 , 0 )-energy equation should give the 
Schrodinger equation in a simple scheme of first quantization. In order to justify such 
an assumption, we must suppose a primordial quantum dynamics of matter contained 
in the universe and the fact that the today observed large scale structure reflects the 
phenomenology of early epochs. In other words, we can say that the Schrodinger Eq. (^) 
describes exactly the behaviour of the universe for any value of the cosmic time and 
the classical Einstein equations (in particular the (0,0)) well reflect the evolution of 
dynamical system when it can be described at a classical level. 

Here we have not considered any questions connected neither to quantum field theory 
nor to quantum gravity. We are simply taking into account the “particle” dynamics of 
objects which are galaxies. The probability of finding them at a certain epoch {i.e. a{t)) 
or a certain red-shift is given by the Schrodinger equation which is nothing else but 
the quantum counterpart of the classical Einstein first order equation. Furthermore, the 
eigensolutions of such an equation, having a probabilistic meaning, can be connected to 
the correlation function which gives the features of the distribution of a certain class of 
objects {e.g. galaxies, clusters or superclusters of galaxies). 

As we have shown, several possibilities exist to get oscillatory behaviours: they de¬ 
pend on the kinds and the number of fluids into the dynamics and, essentially, on the 
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kind of FRW spatial model. If the universe is dust-dominated, with a remnant of cos¬ 
mological constant, the asymptotic behaviour of a Bessel function could easily reproduce 
the BEKS observations (in Budinich 1993, Gegenbauer’s polynomial are used). A good 
match is obtained if '^(a), that is iplz), has 8 oscillations with periodicity of 128/i“^Mpc 
in 2000/i“^Mpc {i.e. in a red-shift range from z = 0 to z = 0.5). The amplitude of such 
oscillations strictly depends on the cosmological densities of fluids involved (see Busarello 
et a/. 1994 for a detailed discussion). 

In conclusion, we can say that a breaking of homogeneity and isotropy at small 
scales (~ 100/i“^Mpc) with oscillating correlations (and anticorrelations) between galax¬ 
ies could be easily implemented considering a sort of cosmological Schrodinger equation. 
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